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Effects of the stellar wind on X-ray spectra of Cygnus X-3 
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ABSTRACT 

We study X-ray spectra of Cyg X-3 from BeppoSAX, taking into account absorption and 
emission in the strong stellar wind of its companion. We find the intrinsic X-ray spectra are 
well modelled by disc blackbody emission, its upscattering by hot electrons with a hybrid 
distribution, and by Compton reflection. These spectra are strongly modified by absorption 
and reprocessing in the stellar wind, which we model using the photoionization code cloudy. 
The form of the observed spectra implies the wind is composed of two phases. A hot tenuous 
plasma containing most of the wind mass is required to account for the observed features of 
very strongly ionized Fe. Small dense cool clumps filling < 0.01 of the volume are required to 
absorb the soft X-ray excess, which is emitted by the hot phase but not present in the data. The 
total mass-loss rate is found to be (0.6-1.6) x IO^^'Mq yr"'. We also discuss the feasibility of 
the continuum model dominated by Compton reflection, which we find to best describe our 
data. The intrinsic luminosities of our models suggest that the compact object is a black hole. 

Key words: accretion, accretion discs - binaries: general - stars: individual: Cyg X-3 - X- 
rays: binaries - X-rays: individual: Cyg X-3 - X-rays: stars. 



1 INTRODUCTION 

Cyg X-3 is a high-mass X-ray binary (XRB) system with a short 
orbital period, P = 4.8 h. It is located at a distance of d ~ 9 kpc 
in the Galactic plane (Dickey 1983, assuming the 8 kpc distance to 
the Galactic Centre; Predehl et al. 2000). Due to strong interstellar 
absorption, its optical counterpart is not observable, though infrared 
observations indicate it is a Wolf-Rayet (WR) star (van Kerkwijk 
et al. 1996). In spite of its discovery already in 1966 (Giacconi et 
al. 1967), it remains poorly understood. In particular, it remains 
uncertain whether its compact object is a black hole or a neutron 
star. Currently, there are two other known X-ray binaries containing 
WR stars, IC 10 X-1 (Prestwich et al. 2007) and NGC 300 X-1 
(Carpano et al. 2007a, b). In both cases, there is dynamical evidence 
that the compact object is a black hole. 

Cyg X-3 is the brightest radio source among X-ray binaries 
(Mccollough et al. 1999). It shows very strong radio outbursts and 
resolved jets (Marti, Paredes & Peracaula 2000; Mioduszewski et 
al. 2001). In the X-rays, it exhibits a wide range of variability pat- 
terns. In particular, transitions between the hard and soft spectral 
state occur on the timescale of months to years. 

The understanding of the radiative processes underlying the 
X-ray spectra of Cyg X-3 remains rather rudimentary. Zdziarski & 
Gierlinski (2004) have shown an overall similarity of the spectral 
states of Cyg X-3 to canonical states of black hole binaries. How- 
ever, details, in particular the energy of the cutoff in the hard state. 
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appear to differ significantly. Many different models have been fit- 
ted to spectra of Cyg X-3 by various authors, but each of them 
appeared to fit the data well, which precluded determination of the 
correct one. The models assumed different system components, ge- 
ometry, radiative processes and the absorbing medium. They in- 
cluded different combinations of absorbers, Gaussian lines, black- 
body, power-law with or without a cutoff, and bremsstrahlung, and 
different model combinations were proposed for the hard and soft 
states, see, e.g.. White & Hoh (1982), Nakamura et al. (1993), Ra- 
jeev et al. (1994). Later, the power-law models were replaced by 
more physical models of Comptonization, see Vilhu et al. (2003), 
Hjalmarsdotter et al. (2004) and Hjalmarsdotter et al. (2008, here- 
after HOB). The interpretation of the intrinsic unabsorbed spectra, 
remains, however, not clear. In particular, H08 show that INTE- 
GRAL data of Cyg X-3 in the hard state could be well modelled by 
three models with very different shapes of the intrinsic continuum, 
yet the model most promising from a physical point of view yielded 
the worst statistical fit. 



A key issue in the search for understanding the spectra of Cyg 
X-3 appears to be the nature of the complex intrinsic absorption, 
most likely caused by the wind from the companion star. In this 
paper, we study the effects of the WR stellar wind on the X-ray 
spectra. In Section|2l we present the data used for our study. Then, 
we discuss our theoretical models in Section [3] In Section |4] we 
describe the complex interactions between the X-rays and the wind, 
and present results of our fits to the data. Sections|5]and[6]contain 
a discussion of our results and our conclusions, respectively. 
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Table 1. The log of the analysed BeppoSAX observations. 



ObsID 


Instrument 


Start date 


Exposure 


Gount rate 






Start time [UT] 


[ks] 


[s-i] 






Soft State 






209530012 


LEGS 


2000 April 9 


11.8 


23.08 ± 0.04 




MECS 


21:02:10.5 


81.0 


14.65 ±0.01 




HPGSPC 




36.0 


12.10 ±0.08 




PDS 


Hard State 


18.0 


2.26 ± 0.04 


209530015 


LEGS 


2000 June 20 


17.5 


4.16 ±0.02 




MECS 


11:05:47.5 


68.9 


5.77 ±0.10 




PDS 




14.0 


28.20 ± 0.07 
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Figure 1. The background-subtracted MECS2 1.65-10.5 keV light curves. 
The orbital modulation atP=: l/5dis clearly seen. 



2 THE DATA 

The BeppoSAX satellite carried on board four co-aligned instru- 
ments, the Low-Energy Concentrator Spectrometer (LEGS; Par- 
mar et al. 1997), the Medium-Energy Concentrator Spectrome- 
ters (MECS; Boella et al. 1997), the High-Pressure Gas Scintil- 
lation Proportional Counter (HPGSPC; Manzo et al. 1997), and the 
Phoswich Detector System (PDS; Frontera et al. 1997). Together, 
they covered the nominal energy range of 0. 1-200 keV, 

We choose here two out of the nine Cyg X-3 observations car- 
ried out by BeppoSAX, one in the soft state (hereafter SS), and one 
in the hard state (hereafter HS). The log of the observations is given 
in Table [T] The LEGS, MEGS and PDS event files were obtained 
from the on-line public archive of the BeppoSAX Science Data Cen- 
tre (hereafter SDC), whereas the HPGSPC files for the SS observa- 
tion have been generated by the BeppoSAX SDC for us. The MECS 
data are available from the units MECS2 and MECS3, which we 
add together. There are no HPGSPC data for the HS observation. 
Following the recommendations of the SDC, we have used the en- 
ergy bands for the instruments as follows: the LEGS, 0.12^.0 keV, 
the MECS, 1.65-10.0 keV, the HPGSPC, 7.0-34 keV, and the PDS, 
15-220 keV. 

The LEGS and MECS data have been selected from circu- 
lar regions centred on the source with an 8' radius. We have used 
blank field observations to obtain the MECS background. Because 
Cyg X-3 is located at a low Galactic latitude of -)-0.8°, and in a 
crowded region, we have used the semi-annuli method (Parmar et 
al. 1999) to obtain the LEGS background. The background for the 
two remaining, high-energy and non-imagining, instruments come 
from off-source observations. The MECS2 light curves are shown 
in Fig.[T] We clearly see strong orbital modulation, characteristic to 
this source. 

The spectra (which are averaged over about 5 orbital cycles) 
have been grouped to satisfy the minimum number of counts in 
each channel as required by the;^^ statistic. The LEGS and MECS 
channels have been grouped up to the minimum of 40 counts per 
bin. The HPGSPC data are grouped by the instrument software ac- 
cording to the signal-to-noise ratio. The PDS spectra have been 
grouped into 18 bins with different minimum number of counts ac- 
cording to the SDC recommendations. 

We then fit simultaneously all of the available instruments. We 
have allowed free normalization of the LEGS and HPGSPC data 
with respect to the MECS data, and the fitted values were about 
0.9-1 for the LEGS, and 0.97-1.05 for the HPGSPC, which is in 
good agreement with the recommendations of the SDC. For the 
PDS data, we assumed the fixed normalization with respect to the 
MECS of 0.86. 




Figure 2. A schematic representation of the geometry of Cyg X-3. The 
spectroscopic inferior conjunction corresponds to the orbital phase of <p = 
0.5. The trajectory of the compact object is plotted in the rest frame of 
the primary. The shaded area represents a spherically symmetric, homoge- 
neous, stellar wind, with the darkness increasing with the density, dependent 
on the distance from the centre of the primary, r. 



3 THE MODEL 

Spectra observed from Cyg X-3 are determined by the form of the 
intrinsic spectra and their modification by the intervening medium. 
The latter consists of absorption in the interstellar medium and ab- 
sorption and re-emission in the local medium, which consists pri- 
marily of the stellar wind. The intrinsic spectrum and the physi- 
cal state of the stellar wind are coupled through complex plasma 
processes. Namely, as the X-ray photons travel through the wind, 
they heat and ionize it, which affects both the wind absorption and 
emission. Thus, when we try to model the observed spectra, we en- 
counter a complication that the observed spectrum depends on both 
the intrinsic spectrum and the distribution of the wind density and 
its physical state along the line of sight, which, in turn, depends 
on the shape and flux of the intrinsic spectrum. Also, the wind ac- 
celeration depends on the irradiating spectrum and the ionization 
state. Below, we deal with some of those issues using the broad- 
band BeppoSAX data. We use detailed physical description of both 
the intrinsic continua and the wind absorption. 



3.1 The binary 

A schematic picture of the binary is shown in Fig.|2] In calculating 
the distribution of the stellar wind, we assume for simplicity that 
the primary is at rest. We also assume that the orbit is circular. The 
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Figure 3. Comparison of the helium WR stellar radii with the Roche-lobe 
radii. The crosses give the model radii of Langer (1989), and the solid curve 
gives the approximation given by equation (20) in that paper. The dotted and 
dashed curves give the Roche lobe radii in Cyg X-3 calculated for the mass 
of the compact object of 1.4Mo and 2OM0, respectively. 
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Figure 4. Wind velocity laws for different values of /Ji and /J2. The solid 
curve represents a fast accelerating wind with /?i = 0, /32 = 1 , the dashed 
curve with j8i = 0, /?2 = 3 is for a slowly accelerating wind, and the dotted 
curve with ji\ = 0.35, jSt = 8 corresponds to winds of some WR stars 
(Springmann 1994). 



mass and the stellar radius of the primary are denoted as Mj, and 
R^,, respectively. The mass of the compact object (either a neutron 
star or a black hole) is denoted as Mx- 

The orbital phase = corresponds to spectroscopic superior 
conjunction (i.e., the compact object being behind the star), and <p is 
within the 0-1 interval. The distance between the centres of mass, 
a, follows from the Kepler law, 

3 _ G(Mx + MJ 



471^ 



3.1 X 10' 



M 



1/3 



cm ^ 4.47 



M 



1/3 



R 



(1) 



(2) 



(30Mo/ \30Mo 

where G is the gravitational constant and M = Mx + M^. 

The primary is a massive WR star, which emits a strong radia- 
tively driven stellar wind. We compare the stellar radii, as given by 
Langer (1989), with the WR Roche lobe radii, calculated follow- 
ing Eggleton (1983) for the separation given by equation l|2j. The 
radii given by Langer (1989) correspond to the star itself, without 
including the wind. Since the wind is, in general, optically thick, 
those radii are smaller than the photospheric radii, i.e., those cor- 
responding to the optical depth of t = 2/3. In Fig. [3] we plot the 
stellar radii corresponding to WR stars with the abundances domi- 
nated by He. As we see, the star always fits well within the Roche 
lobe, regardless of the compact object being a neutron star or a 
black hole. 



3.2 The wind 

We assume the wind to be smooth, spherically symmetric, non- 
rotating, stationary and radiatively driven. We follow Springmann 
(1994) and parameterize the wind velocity as follows. 



v(r) = v^+ (v„ - vo) V/3ix + (l -/3i)^^ 



(3) 



where x = I —R^jr, r is the distance from the centre of the primary, 
v„ is the wind terminal velocity, fi\ and fii are constants which 
parametrize the wind acceleration rate, and vq is the initial wind 
velocity, which we set to a value of 20 km s"', approximately corre- 
sponding to the atmosphere sound speed. According to Springmann 
(1994), equation l[3]l better describes winds of some WR stars than 
the usual wind velocity law of Vo<,(l - Rt,lr^ . 



We consider three different parameterizations of the velocity 
law, also plotted in Fig.|4] 

(i) y6i = 0, j82 = 1 - a very fast accelerating wind; 

(ii) j8i = 0, j82 = 3 - a wind with a slow initial acceleration and 
a fast one at large distances; 

(iii) ySi = 0.35, /32 = 8 - an observational parameterization of 
some WR stellar winds, steep at the base, flat in the middle and 
again steeply accelerating at large distances (Springmann 1994). 

The ion number density (averaged over possible clumping of 
the wind, see in Section [T!4l below) at r is related to v(r) and the 
mass-loss rate, M, via the continuity equation. 



{n(r)) 



M 



(4) 



4m7?H/'o''^v(r)' 

where is the hydrogen mass and yUo is the mean ion molecular 
weight. The wind density along the line of sight is then determined 
by /io, Rt,, M, M, 0, and the inclination, i. 

The elemental abundances of WR stars are significantly dif- 
ferent than the solar. Hydrogen, if present at all, is heavily depleted 
in favour of helium and nitrogen. While the atmospheres of early- 
type nitrogen sequence (WNE) WR stars are hydrogen free, the 
late-type (WNL) ones contain some residual H. Van Kerkwijk et 
al. (1996) classified Cyg X-3 as the WN 4-8 type (i.e., interme- 
diate between those two types) on the basis of the 1 and K-band 
emission lines, but did not find any sign of hydrogen in the infrared 
spectra. Here, we perform and compare calculations for two diff'er- 
ent chemical compositions characteristic for WNE and WNL stars. 
The WR star abundances of H, He, C, N and Fe come from grid 
models of Hamann & Grafener (2004), and are listed in Table [2] 
The remaining abundances were set to the solar values in the sense 
that the mass fractions are set to the solar ones while their relative 
abundances were scaled to the actual abundance of H. 

For a given set of most abundant elements, the H average num- 
ber density is 

{n(r)) 



{nH(r)} 



I + "Hc/H + WC/H + MN/H + "O/H + "Nc/H + "Fc/H ' 



(5) 



where «Hi;/H> "c/h> "n/h> "o/h. WNe/H> "Fe/H represent the abundances 
of the respective elements relative to hydrogen. We set the H mass 
fraction for a WNE star to a small value instead of 0, which would 
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Table 2. The model parameters and the mass fractions and abundances relative to H for the most abundant elements (from Hamann & Grafener 2004) for 
early-type, WNE, and late-type, WNL, WR stars. The H mass fraction in WNE stars is close to zero, here we set it to an arbitrarily chosen small value. 



parameter 


WNE 




WNL 


/J 


1.33 




1.00 






4.05 




2.52 




Vco [cm s~^] 


1.6 X 10* 


1.3 xlO** 




1.78 




1.48 




A>Hc 


1.89 




1.58 




element 


mass fraction 


'^element/H 


mass fraction 


'^element/H 


H 


1.0 X 10-5 


1 


0.2 


1 


He 


0.98 


2.46x10'' 


0.78 


0.98 


C 


1.0 X 10-"* 


0.84 


1.0 X lO-'' 


4.2 xlO-5 


N 


1.5 X 10"- 


1.08 xlO^ 


1.5 X 10-^ 


5.4x10-3 


O 


5.5 X 10-3 


34.65 


5.5 X 10-3 


1.73 X 10-3 


Ne 


1.42 X 10-3 


7.07 


1.42 X 10-3 


3.54 X 10-** 


Fe 


1.4 X 10-3 


2.53 


1.4 X 10-3 


1.26x10-^* 



otherwise prevent us from defining the abundances with respect to 
H. The abundances used by us are listed in Table |2] WNE stars 
tend to have faster winds than WNL ones (Hamann, Koesterke & 
Wessolowski 1995). Therefore, we choose the velocity range of 
= 1300-1600 km s-' (van Kerkwijk et al. 1996), and attribute 
the lower limit to WNL stars and upper limit to WNE stars. The 
wind is likely to be not smooth, but clumped, which effect we in- 
clude in the way described in Section [34l below. 

3.3 The intrinsic X-ray spectra 

In spite of decades of research, the form of the intrinsic spectra 
of Cyg X-3 remains highly uncertain. The reason for it is the very 
strong X-ray absorption in this source. As mentioned above, the 
state of the absorbing medium is coupled to the form of the ir- 
radiating spectrum. This introduces strong coupling between var- 
ious fit parameters. The X-rays are further absorbed in the inter- 
stellar medium and reach the observer with almost no photons be- 
low 1 keV. Only the high energy part of the spectrum remains rela- 
tively unabsorbed. In order to efficiently deal with this problem, we 
start with a simple phenomenological description of the absorption 
while fitting the parameters of the intrinsic spectrum. 

The X-ray radiation, most likely, originates in the vicinity of 
the compact object, and it is produced via accretion. The accretion 
disc emit soft blackbody photons. Those photons are Comptonized 
in a hot plasma in the vicinity of the disc. The plasma contains en- 
ergetic electrons, which distribution may be hybrid, i.e., containing 
both a Maxwellian part and a nonthermal tail beyond it. Similarly to 
Vilhu et al. (2003), Hjalmarsdotter et al. (2004) and H08, we model 
the continuum emission of such a plasma with the eqpair model 
(Coppi 1999; see also Gierlinski et al. 1999), as implemented in 
the xspec spectral fitting software (Amaud 1996). The model cal- 
culates the spectrum from a spherical plasma. The distributions of 
photons and electrons are assumed to be isotropic and homoge- 
neous. Some major model parameters are defined via the compact- 
ness. 



where L is a luminosity, R is the characteristic size of the plasma, 
and (Tx is the Thomson cross section. In particular, we define the 
compactness, corresponding to the soft photons irradiating the 



hot plasma with the luminosity, L,. The soft photon spectrum is 
taken here to be a disc blackbody (diskpn of Gierlinski et al. 
1999), which is characterized by the maximum blackbody temper- 
ature, Tbb- The normalization, Nc, of the eqpair model equals that 
of the input disc emission, diskpn (Gierlinski et al. 1999). 

The electron distribution contains a Maxwellian part, with the 
optical depth corresponding to ionization electrons, Tp. The thermal 
electrons are heated, i.e., supplied with the power, Lth- If the plasma 
contains also a nonthermal part, it is assumed that the electrons are 
accelerated at a power-law rate oc /-"^'"J, where y is the Lorentz 
factor and Finj is the acceleration power-law index. The electrons 
are accelerated out of the thermal distribution, but in steady state 
the same number of electrons returns to it after losing their en- 
ergy. The total power in the acceleration is L„,i„ and the total power 
supplied to the plasma is Lh = Lfa + Lmh- We then define three fur- 
ther compactness parameters, ^th, imh, and 4, according to equation 
l|6j. The physical processes taken into account are nonthermal elec- 
tron energy losses via Comptonization and Coulomb losses with 
the thermal plasma, cooling and heating of the thermal plasma by 
Comptonization, and production of e* pairs. A process also taken 
into account, but usually negligible in luminous accretion flows, 
is bremsstrahlung. Then, the steady-state electron distribution, in- 
cluding the temperature of the thermal part, T^, is solved for self- 
consistently, and the corresponding emitted photon spectrum is cal- 
culated. The total steady-state optical depth may be higher than Tp 
due to e* pair production. The overall hardness of the spectrum is 
determined by the ratio, 4/4, i-e., the amplification factor of the 
Comptonization. 

A fraction of the X-rays is scattered back to the optically- 
thick medium emitting blackbody photons (presumably the ac- 
cretion disc), and it is reflected with an energy-dependent albedo 
(Magdziarz & Zdziarski 1995). The relative strength of reflection 
is parametrized by Q./2n, where SI is the effective solid angle sub- 
tended by the reflector. If the primary emission is isotropic and nei- 
ther the primary emission nor the reflected one is obscured or atten- 
uated (apart from any external absorption, e.g., in the wind and the 
interstellar medium), SI is also the actual solid angle subtended by 
the reflector. However, the intrinsic X-ray source may be obscured, 
e.g., by an inner torus-like structure, whereas the reflection from 
the torus back wall can be seen directly. Then, n/2;r is a measure 
of the degree of obscuration of the central source. The reflecting 
surface is allowed to be ionized, with the ionization parameter de- 
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fined as ^ = iion/«refl^refl. where Li„„ is the 5 eV-20 keV luminosity 
and fjieti, r,cfl is the density and the temperature, respectively, of the 
reflector located at a distance, R^^f], from the illuminating source. 
We use here the ionization model by Done et al. (1992), which al- 
though is not applicable to highly ionized Fe (Ballantyne, Ross & 
Fabian 2001), is sufficiently accurate at low ionization, such as that 
found for our data below. We hereafter denote the total intrinsic 
continuum, which includes the disc blackbody emission and both 
Compton upscattering and Compton reflection, as Fq ■ 

The primary and reflected emission is then absorbed in the 
intervening media. One of the absorbers is the interstellar medium, 
with the line-of-sight column density, A'ho- It has been estimated, in 
units of 10^" cm-^ as 1.7, 1.4 and 1.42 by Chu & Bieging (1973), 
Lauque, Lequeux & Nguyen-Quang-Rieu (1972), and Dickey & 
Lockman (1990), respectively. Here, we assume A'ho = 1-5 x 
10^^ cm"^. We assume that absorbing medium is neutral, model it 
with the wabs model (Morrison and McCammon 1983), and denote 
the resulting energy-dependent attenuation factor as Aq(E), where 
E is the photon energy. 

In addition, there is a strong local absorption, apparently in 
the stellar wind, as none of our data can be satisfactorily fitted with 
interstellar absorption only. However, as discussed above, the in- 
trinsic continuum and the absorbing wind parameters are strongly 
mutually dependent. Thus, we start with a simple model of intrinsic 
absorption modifying the eqpair continuum. Namely, we assume 
a product of complete and partial covering, in both cases assum- 
ing neutral absorbing medium. The column densities are A'n.tot and 
A^H.part, respectively, and the second absorber covers a fraction, C, 
of the source, 

ApJE) = [ce-'^<'^>"H,p.- + 1 _ c] e-"^"^*""-'-', (7) 

where (t(E) is the absorption cross section. This phenomenological 
model can fit the data relatively well, and was used for Cyg X- 
3 spectra by Vilhu et al. (2003), Hjalmarsdotter et al. (2004) and 
H08. 

Furthermore, the data show a relatively strong line related to 
the Fe K complex. We model it as a single Gaussian with the cen- 
troid energy, Ep^, the width, crpc and the normalization, /pc. The 
iron lines originate both from reflection and the local absorber. We 
denote this component as Fpc(E). 

The total model has then the form, 

F{E) = Ao(£)Ap„-(£) [FciE) + Fp,{E)] . (8) 

Its main free parameters are the Comptonization amplification fac- 
tor, ^h/4> the non-thermal fraction, ^nih/4> Tp, kTf,f„ the reflec- 
tion strength, Sl/2n, T^^j, A'n.iot and A'H.pan and C. Then, we as- 
sume fixed values of 4 = 100, r,cfl = lO'K (important for re- 
flection only), the characteristic size of the plasma (important only 
for the Coulomb and bremsstrahlung processes), R = lO' cm, the 
minimum and maximum Lorentz factors the accelerated electrons. 
Train = 1-3, 7raax = 1000, respectively, and the inclination angle, 
i = 60° (important for the disc emission and reflection). We neglect 
relativistic broadening of the reflected spectrum and the Fe K line. 

The relative iron abundance, Ape, and the relative abundance 
of elements heavier than He, A>hc, used by eqpair, assume the 
solar elemental abundances of Anders & Ebihara (1982). They are 
important for reflection, and determine the relative importance of 
bound-free absorption (determined in the X-ray regime mostly by 
elements heavier then He) to scattering (determined by the electron 
abundance). We rescale them with respect to the actual electron 
abundance in the wind with the WR-star abundances, see Table |2] 
Namely, 



(«Fi;/«i;)wR . (n>Hc/«c)wR ,„s 

Ape = -— , A>He = -— , (9) 

("Fe/nc)o ('l>He/'1c)G 

where the WR and solar subscripts refer to the WR-star abundances 
used by us and to those of Anders & Ebihara (1982), respectively, 
and, for each, the ratios to are given by, 

"Fe ^ "Fe/H (Mo/f^-^T^ ^jq^ 

«c 1 + «Hc/H + MC/H + "N/H + no/H + "Ne/H + "Fe/H ' 
">He _ ("C/H + "N/H + "O/H + "Ne/H + "Fc/h) (Mo/t^ - 1)~' ^^^^ 
«c 1 + «He/H + «C/H + «N/H + «0/H + «Nc/H + «Fe/H 

where fj. is the mean molecular weight. 
3.4 Radiative transfer in tlie wind 

The X-ray spectrum emitted by the hot plasma is then transferred 
in the stellar wind. This radiative transfer is modelled using the 
cloudy photoionization code (v. 6.02, see Ferland et al. 1998). The 
main input parameters of cloudy are the shape of the irradiating 
continuum, the density profile along the line of sight, and the wind 
elemental abundances. The code also allows to account for the ef- 
fect of possible wind clumping, with clouds filling a fraction, /, of 
the volume, and with empty region between the clumps. 

For a given M, the density, n, of the matter inside the clumps 
equals {n)/f, i.e., it is 1 // times higher than that of a smooth wind, 
{n}. Although neither (n) nor the column density change, the state 
of the plasma does. In photoionized plasmas, the volume-averaged 
ionization rate is oc ((«)//)/ = (n). On the other hand, the recom- 
bination rate is proportional to n^, and its volume-averaged rate is 
ff'f = {nf' I f, i.e., it increases. Thus, for a fixed radiation 
field, increasing the density decreases the average effective charge 
of a given element, i.e., the plasma becomes less ionized. 

In our calculations, we also consider a two-phase wind, in 
which the region between the clumps is not empty. We approxi- 
mate that situation by considering two cloudy models, one with 
a smooth, f = \, optically thin wind, and one with the clumps, 
/ <sc 1. The results of the two calculations are then superimposed 
on the intrinsic spectrum, see Section l42l below. 

We note that in our calculations we neglect the effect of X- 
ray irradiation on the wind velocity and density. Such calculations 
cannot be done within the framework of our model. They would re- 
quire iterative hydrodynamic calculations of the entire wind struc- 
ture, not only the effect of the irradiation on the line of sight. Since 
we neglect the effect of the X-ray radiation force on the wind, we 
consequently do not consider the pressure balance in the acceler- 
ated wind, which would taking into account the pressures of radia- 
tion and each phase of the wind as well as of the wind inertia. We 
also point out that cloudy does not allow for the filling factor with 
an arbitrary dependence on position, and thus it would not allow us 
to achieve a global pressure balance of the two wind phases. 

We use the wind density distribution along the line of sight 
from the X-ray source calculated for a wind spherically symmetric 
with respect to the primary star. This density distribution is then 
used by cloudy to create a cloud which is spherically symmetric 
with respect to the X-ray source, an approximation required by that 
code. Also, all of the radiation scattered from (or reemitted by) one 
side of the wind back to the opposite side is included in calculations 
of the ionization structure and line emission (which is called the 
closed geometry in cloudy). It is also important to note that in the 
case of closed geometry, cloudy does not calculate the effects of 
electron scattering, which may be important for electron column 
densities substantially greater than 10^^^ cm"^. 
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Figure 5. The dependencies oi xl on kT^,^, in tiie HS and SS for the phe- 
nomenological treatment of absorption and the WNE stellar type (see Sec- 
tion |4.U . The number of degrees of freedom is 466 and 521 in the HS and 
SS, respectively. 



For the irradiating continuum, we use the results of the first 
stage of the calculations assuming the phenomenological absorp- 
tion, equations l|7}{8). We create a set of table models in the xspec 
format (both additive and multiplicative) for the HS and SS spectra 
and for a number of the wind parameters, which we then use for the 
final spectral fitting. 

One issue we need to consider here is the effect of the black- 
body photons emitted by the primary on the wind. Those UV pho- 
tons are important for the balance between Compton heating and 
cooling. In the limit of complete ionization, the gas reaches the 
Compton temperature, Tc, 



f EF{E)dE 
4k / F(E)dE' 



(12) 



(e.g., Begelman, McKee & Shields 1983), where the irradiating 
flux, F, is the sum of the X-ray flux. Fx, and the stellar emission. 
Ft . The numerator and denominator in equation l ll2t correspond to 
the heating and cooling, respectively. We note that for two separate 
sources of illumination, the ratio between the two fluxes is a strong 
function of the position in the wind, except for large distances from 
the binary, where, in the optically thin case, it reaches a constant. 
That approximation, with the sum of the X-ray and UV sources as 
ionizing the wind, was appropriate for absorption far away from 
the binary plane considered by Szostek & Zdziarski (2007) for Cyg 
X- 1 . However, it cannot be applied to wind absorption close to the 
binary plane, as it is the case in our present study. 

A luminous WR primary may be a significant contributor to 
the cooling. On the other hand, the X-ray source with Lx ~ 10'* 
erg s"' (Vilhu et al. 2003) dominates the heating except very close 
to the primary. However, the cloudy code does not allow to in- 
clude two independent and spatially resolved illuminating sources. 
Therefore, we consider the X-rays as the only source of the wind 



hard state 




1 5 10 50 100 

Energy [keV] 

Figure 6. The top panel shows the best-fit unfolded spectra (crosses) and 
the best-fit models (solid curves) obtained with the phenomenological treat- 
ment of absorption, see Section l4~Tt . The dashed and dotted curves show 
the unabsorbed spectra in the HS and SS, respectively. The two lower pan- 
els show the corresponding data to model ratios. The black, red, blue and 
green symbols correspond to the LEGS. MECS, HPGSPC, and PDS data, 
respectively. 



illumination. Since we neglect the existence of the stellar cooling, 
the model gives an upper limit on the wind temperature. Thus, it is 
somewhat more ionized, and less absorbing, in our model than in 
reality. Hence, it gives upper limits on the wind density and M for 
a given set of model parameters. In Section |4~8l below, we estimate 
the maximum effect of the wind illumination by the companion on 
our final results, and find it very minor. 



4 RESULTS 

4.1 The phenomenological treatment of absorption 

We first fit the data using the phenomenological treatment of ab- 
sorption as a product of complete and partial covering, equation 
([8}. Still, even this simplified model shows a complex structure, 
with a number of local minima and strong correlations between its 
parameters. In particular, the values of the maximum blackbody 
temperature, Tbb, and of the absorber column densities are strongly 
correlated. The absorption removes almost all photons below 1 keV, 
which results in very weak constraints on the form of the intrinsic 
spectrum in that range. For very similar fits to the data, those in- 
trinsic spectra can have very different form. This uncertainty would 
also affect the ionization structure of the wind (considered in Sec- 
tionl42lbelow). 

Thus, we first perform a detailed search for the global min- 
imum as a function of r(,b- Since the HPGSPC and PDS spectra are 
only weakly affected by the absorption whereas the lower energy 
spectra are strongly affected, we temporarily increase the relative 
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weight of the HPGSPC and PDS data by decreasing their measure- 
ment errors by a factor of 10. Also, given the relatively poor quality 
of the high energy data in the SS, we use for it the constant T^^j = 2. 
Also, the Fe K line complex is relatively weak in the SS, and we 
thus fix CTpj. = 0.2 keV. We also find that changing the abundances 
between the WNE and WNL models (Table |2j has little effect on 
the intrinsic spectra. Thus, we consider here only the WNE model. 

Because the data are averaged over ~5 orbital periods, it is 
unlikely that all discrete spectral features would be fully correctly 
reproduced with a single-phase model. Therefore, we add a 3 per 
cent systematic error to the model. Because of a large number of 
local minima, the default xspec minimization method of modified 
Levenberg-Marquardt based on the curfit routine from Beving- 
ton & Robinson (1992) was found to be inefficient. Thus, we use the 
migrad method based on the CERN minuit code (James & Win- 
kler 2004), which is also implemented in xspec. The latter method 
is more time consuming but was found to give much better results 
than the former. 

We present the resulting dependencies of xl on kTb\, in Fig. 
|5] We see a relatively narrow minimum in the SS at ~1.2 keV. On 
the other hand, the minimum in the HS is much broader, at ~1.5- 
1.9 keV. Hereafter, we assume constant values of kT^b = 1.2 keV 
and 1.7 keV in the SS and HS, respectively. The higher value of 
kTf,], in the HS may be due to the presence of some kind of a soft 
excess in that state, commonly found in black-hole binaries (e.g., 
Gierlihski & Zdziarski 2005), or may be due to the origin of at least 
some of the soft photons from another form of matter, e.g., dense 
blobs in the wind (see Section l4!2l below). The HS and SS unfolded 
spectra, including the unabsorbed, intrinsic, spectra are shown in 
Fig. |6] (after restoring the original HPGSPC and PDS errors). The 
HS and SS spectra studied here are very similar to the spectra 1 
and 5, respectively, obtained from averaging pointing observations 
of Rossi X-Ray Timing Explorer by Szostek & Zdziarski (2004). 

The rather high minimum values of xl in both states point to 
the model being too simplified. Its main problems can be identified 
by looking at the fit residuals in Fig. [6] (with the original HPGSPC 
and PDS errors). Strong residuals at energies < 3 keV are appar- 
ently from neglecting the wind emission lines. Also, the Fe K com- 
plex is not accurately treated by a single Gaussian. Still, the overall 
continuum shapes have been found to be sufficiently accurate to 
calculate the effect of the wind irradiation, see below. 

4.2 Wind absorption and emission 

We use the above best-fit unabsorbed models as the input spectra 
for the cloudy calculations of the radiative transfer in the stellar 
wind. The 13.6-2 x 10^ eV ionizing luminosities (used by cloudy) 
are 8.0 x 10" erg s"' and 2.8 x 10'* erg s"' in the HS and SS, 
respectively. 

Before the actual fitting, we calculate the effect of changing 
the wind parameters on the transmitted spectra. Since the fitted 
spectra are phase-averaged, we use the orbital phase in the mid- 
dle between the superior and inferior conjunctions, = 0.25. This 
choice of has also the advantage that, assuming spherical sym- 
metry, there is no dependence of the wind absorption on the in- 
clination. (The only effect of i is then on the disc blackbody and 
reflection components of the intrinsic continuum.) For these calcu- 
lations, we consider the HS only but the conclusions also apply to 
SS. Then, we assume as the test default parameters the WNE stellar 
type, the velocity law (i) (with /?, = 0, yS^ = 1), R^ = I x 10"cm, 
the total mass of the two components, M = 5OM0, / = 1 and M = 
2.0 X 1O"''M0 yr"'. For these parameters, the column density of 
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Figure 7. Example model spectra transmitted through the stellar wind for 
a range of values of M (colour curves) and the rest of the parameters kept 
at the values given in Section |4.2| The original unabsorbed continuum is 
shown by the black curve, (a) The spectra transmitted only through the stel- 
lar wind, (b) The spectra transmitted through both the wind and the inter- 
stellar medium in the direction of Cyg X-3. 

(mostly fully ionized) He along the line of sight is A'hc = 4.2 x 10^^ 
cm"^. The corresponding Thomson optical depth is r„ 0.6. We 
change one of the default parameters while keeping the rest at the 
above values. 

Fig. Ha) shows the spectra transmitted through the wind for 
M = (1.6-5.0) X 1O"'M0 yr"'. As expected, the increasing M 
results in the increasing absorption. The wind is strongly ionized 
and, except for the two highest values of M, there is almost no ab- 
sorption at energies < 1 keV. On the other hand, there is relatively 
strong wind emission at those energies. Fig. |7lb) shows the same 
spectra but transmitted also through the interstellar medium with 
A'h of Cyg X-3. Now, all the emission at < 1 keV is strongly ab- 
sorbed. We would like to draw the attention to the shape of the 
spectrum for M = 2.4 x IO^^Mq yr"'. The elements which con- 
tribute to absorption below 1 keV are completely ionized and the 
low energy photons are not absorbed. The elements responsible for 
absorption above 1 ke V are only partially ionized and create a deep, 
broad absorption feature. The combined effect of low absorption in 
the wind at low energies and strong interstellar absorption creates 
a spectral shape with a strong soft excess at < 2 keV. 

Fig. [8] shows the effect of changing the other parameters on 
the spectra absorbed by both the wind and the interstellar medium. 
The wind of the WNL-type star absorbs much more than that of the 
WNE type due to the lower wind velocity of the former, Fig.[8lb). 
Absorption decreases with the increasing total mass (determining 
the separation) and the stellar radius (which is due to the scaling 
of the wind velocity with r/Rt,, see Section [J!2t . see Figs.[8la, d). 
Then the wind velocity laws (i)-(iii) yield correspondingly increas- 
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Figure 8. Example model spectra transmitted through both the wind and the interstellar medium. In each of the (a-e) panels, one parameter is varied (identified 
by its colour) while the rest are kept at the initial values, given in Section 132] The black curves coirespond to the initial value of the variable parameter, (t) 
Coordinated variability of / and M. 



ing absorption, Fig.[8fc). In most cases, the effect on the spectra at 
~ 1.5 keV is only minor, as the elements responsible for absorption 
in that range are almost fully ionized in the wind. 

Fig-Hie) shows the case where a coordinated decrease of both 
/ and M leads to the overall form of the spectrum at > 2 keV chang- 
ing only slightly. On the other hand, the reduction of / leads to a 
decline and disappearance of the emission excess at < 2 keV. The 
lower the /, the higher the density (of the blobs), and the weaker 
the ionization, reducing then the soft X-ray emission. Another ef- 
fect we see is the change in the spectral shape in the ~ (6-10) keV 
region. Here, the hotter, more ionized, gas (at a higher /) produces 
stronger Fe emission and absorption than the cooler one at lower /. 



4.3 The minimum and maximum mass-loss rates 

Then, we would like to obtain the lower and upper limits on M, 
which correspond to the model parameters maximizing and min- 
mizing (for a given M), respectively, the absorption. The minimum 
corresponds to the WNL-type star, the (iii) velocity law, and the 
smallest possible separation between the system components. For 
the last requirement, we would like to consider the smallest possi- 
ble masses and the highest stellar radius. We assume Mx = 1.4M0, 
= 5Mq (which is close to the smallest possible WR-star mass, 



and agrees with the estimates of Hanson, Still & Fender 2000), 
which corresponds to = O.6R0 (Fig.[3]l. 

Next, the maximum of M corresponds to the WNE-type star, 
the (i) velocity law, and the largest possible separation between the 
system components. For the last requirement, we would like to con- 
sider the largest possible masses and the lowest stellar radius. We 
assume Mx = 20Mq, M-,, = 50Mq (which is close to the largest 
possible WR-star mass using the standard models of Langer 1989), 
which corresponds to R^, = 2.3Ro. (We note that the assumed black 
hole mass is larger than that allowed by the constraints of Hanson 
et al. 2000). 

For both the minimum and maximum we consider a number of 
values of /. For each /, we create xspec tables for a range of val- 
ues of M. We will denote the resulting wind absorption as A„(£), 
and the corresponding line emission, as F„(£'). The following con- 
siderations apply to both minimum and maximum limit on M, but 
we plot the results only for the maximum case. 



4.4 Homogeneous hot wind 

We first consider the case with no clumps, / = 1. As before, we 
keep the constant Finj = 2 in the SS, and kT^b = 1.7 keV and 1.2 
keV in the HS and SS, respectively. We also add a 3 per cent sys- 
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Figure 9. The best-fit unfolded spectra and residuals fitted including wind absorption and emission (for the case maximizing M). The colours of the residuals 
are the same as in Fig.|6] The HS fitted by the model of equation (13), CD, and {BJ are shown in (a), (c) and (e), respectively. The SS fitted by the model of 
equation )13t . (14), and (15) are shown in (b), (d) and (f), respectively. The emission of the hot wind, / = 1, and the cool one, / « 1, are shown by the cyan 
and blue lines, respectively. 



tematic error to the model (as in Section \4J[ and a 1.5 per cent 
systematic error to each data set. We drop other simplifying as- 
sumptions, in particular absorption is now only in the wind and the 
insterstallar medium. The structure of the model is, 

F{E) = Ao{E) [A„,h(£)Fc(£) + Fw,h(£)] • (13) 

where the subscript 'h' indicates that the wind with / = 1 is hot. 
Below, we also use the subscript 'c' for the cooler gas in the com- 
pressed clumps with / «; 1. In the spherical geometry with the 
central irradiating source used by cloudy, the relative normaliza- 
tion of the emission-line component is unity. However, we use here 
cloudy only as an approximation to the actual stellar wind geom- 
etry. Thus, we allow a free relative normalization of the emission- 



line component, but we find its best fit values to be ~1. We denote 
that factor by Wh, for the cases with / = 1, / 1, respectively. 

We show the resulting fits and their residuals in Figs. |9la,b). 
We see the model provides rather bad fits, mostly due to its predic- 
tion of an excess emission at low energies. As explained in Section 
14.21 this excess is created by low absorption by the wind below ~2 
keV, see Fig.|7] The wind emission gives also some contribution, 
see the cyan lines in Figs.[9la,b), but we can see it is rather minor 
(in comparison to the total model). We have tested this by fitting a 
model with the wind emission removed, and obtained rather similar 
soft excess in the model compared to the data at < 2 keV. Also, the 
SS model strongly overestimates the depth of the Fe xxv and xxvi 
edges at ^9 keV (see, e.g.. Turner et al. 1992 for the Fe K edge 
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Figure 10. The best-fit absorbed HS model (for the case maximizing M) 
including two-phase wind and an Fe Ko- line from Compton reflection (top 
solid curve). The wind emission is plotted in green below the top curve. The 
fit residuals are plotted in the same colours as in Fig.|6] 

and line energies). Thus, the data require the presence of additional 
denser material with lower ionisation. 

4.5 Clumpy cool wind 

In order to deal with the < 2 keV residuals, we can reduce the filling 
factor, /, see Fig.[8le). From our fitting, we find that the maximum 
values of / that remove those residuals are approximately / = 0.01 
and 0.005 in the HS and SS, respectively. The model is, 

F(E) = Ao(£) [A„,c(£)Fc(£) + F„,,(£)] . (14) 

We show the resulting fits together with their residuals in Figs. 
[9lc, d). We see that this model provides a better fits than that with 
/ = 1, and ;t'2 ^ j (479 d.o.f.) and 1.97 (524 d.o.f.) for the 
HS and SS, respectively (for the case maximizing M). Still, there 
are significant residua at various energies, in particular there are 
signs of highly ionized Fe in the data, which are not included in 
the model with a low /, in which the wind clumps are relatively 
weakly ionized. This indicates the presence of a hot phase. 

4.6 Two phase wind 

Indeed, it is highly unlikely that there is only vacuum in the inter- 
clump wind medium. Therefore, we now model the wind as con- 
taining both the cool clumps and a smooth hot gas between them. 
Because the two phases coexist together in the same space, as dis- 
cussed above, we account for it as a combination of two separate 
cloudy models calculated independently for the same system pa- 
rameters and the illuminating source but for different M and /. This 
approximation is valid as long as the hot phase is optically thin, i.e., 
it does not substantially absorb soft X-rays but still emits in lines. 
The hot, very strongly ionized, phase may still cause some absorp- 
tion around ~9 keV. Each of the phases is now characterized by its 
own accretion rate. The model is now, 

F{E) = Ao{E) [A„.,(£)A„.h(£)Fc(£) + F„.,(£) + F„.h(£)] . (15) 

We show the resulting fits together with their residuals in Figs. 
[9le,f)- The overall fit quality is now very good, with;^'^ = 1.03, 0.99 
(468 d.o.f.) in the HS at the minimum and maximum M, and^J = 
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Figure 11. The unabsorbed spectra (solid curves) of the emission before 
the wind reprocessing shown together with their components for (a) the HS 
and (b) the SS. The unscattered disc blackbody, Compton upscattering by 
thermal and by nonthermal electrons, and the Compton reflection together 
with the Fe Ka line are shown by magenta dot-dashes, red short dashes, 
green long dashes, and blue dots, respectively. The cyan solid curves show 
the intrinsic spectrum without the reflection component. 

1.20, 1.29 (522 d.o.f.) in the SS at the miniumum and maximum 
M, respectively. The accretion rate in the hot phase is always larger 
than that in the clumps. 

We note the presence of an absorption feature at ~4.4 keV 
in the HS MFCS data, but not in either of the LEGS data. This 
energy corresponds to an Xe absorption line and it is most likely 
of instrumental origin (SDC, private communication). Adding an 
absorption line at that energy (which we keep for the following fit) 
improves the fit to;^'^ = 0.99, 0.96 (465 d.o.f.) at the minimum and 
maximum M, respectively. 

The two-phase model fits the SS data slightly worse than the 
HS ones. The SS xl 's relatively large mainly due to the residuals 
around 1.6 keV. By including a narrow line at that energy into the 
model, we were able to considerably improve the fit quality. The 
minimum and maximum M xl become 0.95 and 1.06, respectively 
(519 d.o.f.). We suspect that that line could be produced by Fe xxiii. 
However, since we have no certainty about its origin, we leave it out 
of our final model given below. 

4.7 Compton reflection, the Fe K hne and the final models 

The HS model is extremely reflection dominated, with n/2;r reach- 
ing 50, which can be achieved if the X-ray source is obscured and 
we only see the light reflected from some optically-thick medium. 
(Note that in our models, the wind is illuminated by the spectrum 
fitted in Section |4n which is also reflection-dominated.) Compton 
reflection should also give rise to an Fe Ka line, which is neither 
included in the eqpair model nor its presence is apparent in the 
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residuals. Its equivalent width with respect to the reflection spec- 
trum is ~1 keV (Zycki & Czerny 1994) for hard incident spectra, 
with the photon index < 2. Therefore, we add a Gaussian line to 
the model and fit again the data but imposing the condition that the 
equivalent width, EW = 1 keV, with respect to the reflection com- 
ponent. This lowers the best-fit value of Q./2n. This model has the 
form, 

F(E) = (16) 
Ao(£) jA„.,(£)A„.h(£) [Fc{E) + Fpe(£)] + F„.,(£) + F^,JE)] . 

Adding the line improves the fit significantly for the HS data, the 
F-test probability for the improvement being by chance (Bevington 
& Robinson 1992) is < 3 x 10"^ (If the Fe line is allowed to have a 
free normalization, it changes only slightly with respect to the fixed 
EW = 1 keV.) Note that the best-fit line energy is ~6.7 keV, and it is 
moderately broad, o-pc ~ 0.3 keV. The fit parameters are given Ta- 
ble[3] and Fig.[TO]shows unfolded spectrum with the residuals. Fig. 
1 1 1 f a) shows decomposition into the model components (following 
the method described in Hannikainen et al. 2005). 

The same procedure in case of SS was not successful. It was 
not possible to introduce an additional Fe line to the model (al- 
though its predicted equivalent width is ~3 times lower because of 
the softness of the spectrum above 7 keV). Thus, we keep the model 
as described in Section |4~6l and plotted in Fig. Of)- The best-fit pa- 
rameters are given in Table[3] Fig. II li b) shows decomposition into 
the model components. In all of the models, the ionization param- 
eter, ^, was found to be <k 1 , consistent with neutral reflection. 

Fig. [12] shows the resulting wind temperature profiles as a 
function of the distance from the X-ray source. We see the wind is 
strongly heated by the X-ray source, achieving temperatures much 
higher than that of the stellar radiation. Close to the X-ray source, 
the wind is at the Compton temperature, equation JI2b . i.e., very hot 
and almost completely ionized. Far away, both the wind and photon 
density decrease proportionally to the square of the distance. Thus, 
each phase reaches a constant ionization parameter, ^„ = L/(nr^). 
Since the wind is then optically thin, its state is determined only 
by (Tarter, Tucker & Salpeter 1969), and then the temperature 
of each phase approaches constant, as seen in Fig. [T2] The wind 
Thomson optical depth along the line of sight is 0.31 and 0.43 in 
the HS and SS, respectively. Because cloudy does not take elec- 
tron scattering into account in radiative transfer calculations, the 
resulting intrinsic spectra corrected for scattering will have higher 
fluxes (see Yaqoob 1997) than those given in Table [3] While the 
effect is relatively small in our case (~10 per cent), it can be signif- 
icant around the superior conjunction. 



4.8 Cooling by the primary 

As the final step, we now estimate a possible change of the wind 
parameters caused by the WR stellar photons, see Section [34l Sin- 
gle WNE/WNL WR stars have ^ (1-30) x 10'** erg s"' and the 
effective temperatures of T^g (3-9)x 10* K (Hamann et al. 1995). 
To estimate the maximum cooling eff^ect we add the blackbody of 
the brightest and coolest star (L* = 3 x 10'' erg s"', Teff = 3 x 10" 
K) to the intrinsic continuum assuming they both originate from 
the position of the compact object. This maximizes the effect of the 
blackbody radiation on the plasma in the line of sight. We then run 
the radiative transfer calculations to create emission and absorption 
models for the SS case with the maximum M. 

We find that the inclusion of such strong cooling source in- 
fluences only slightly the wind parameters. The M of the hot com- 




Distance from X-ray source [cm] 



Figure 12. The wind temperature as a function of the distance from the 
X-ray source. The solid and dotted curves conespond to the homogeneous 
wind and the clumps, respectively, in the hard state. The dashed and dot- 
dashed curves correspond to the homogeneous wind and the clumps, re- 
spectively, in the soft state. 

ponent decreases by only ~1 1 per cent and becomes 9.3 x lO^^M© 
yr"' . The cool clumpy phase is not affected by the addition of black- 
body photons and its M remains within its fit uncertainties given in 
Table [3] Thus, not including the stellar blackbody in our calcula- 
tions has a negligible effect on our results. 

5 DISCUSSION 

5.1 The mass-loss rate 

An important result of this study is a new estimate of the mass loss 
rate independent of those of most of other Cyg X-3 studies. The 
M for either HS or SS has been found to be in the range of (0.6- 
1.6) X 10"' Mq yr"'. If the terminal wind velocity is different, v'^, 
from that assumed by us, Vo<, = (1.3-1.6) x 10** cm s"', our values 
of M should be rescaled by ~ v'^lv^. Our values agree with the 
usual range of values of M for WR stars of (0.5-15) x IO^'Mq yr"' 
(Hamann et al. 1995). 

On the other hand, we disagree with Mitra (1996), who, sim- 
ilarly to us, estimated absorption in the WR star wind and claimed 
that for M > 1 x lO^'M© yr"' the wind will be opaque to X- 
rays. He based that limit on some approximate considerations of 
the wind opacity, and we are unable to fully explain the discrep- 
ancy with our detailed quantitative results using the photoioniza- 
tion code cloudy. A factor contributing (but not fully explaining) 
to this difference is his use of lower values for both and the 
separation between the stars. This leads to an increase of the den- 
sity at the location of the compact object by an order of magnitude 
with respect to our calcuations. He did use the cloudy model in a 
single case with a rather high M = 6.3 x 10"' yr"', for which 
his finding of the wind being completely opaque agrees with our 
results. 

Kitamoto et al. (1995) attributed the orbital period change to 
the mass loss from the binary. They found, from this requirement, 
M = 5.8 X 10"'(M/Mq)Mq yr"' . For consistency with our results, 
M ^ (10-30)Mq is required. This includes a large range of possi- 
ble masses of WR stars as well as allows for the presence of either 
a black hole or a neutron star. 

Then, Waltman et al. (1996) found M< 1.0 x IQ-'M© yr"' 
from the observed time delays of post-outburst jet peaks at lower 
frequencies with respect to higher frequencies. Such delays set lim- 
its on the opacity of the stellar wind. Their limit overlaps with our 
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Table 3. The fit parameters of the best-fit models, equation U5K for the cases minimizing and maximizing M. An additional Fe Ka line from Compton reflection 
with the EW = 1 keV with respect to the reflection spectrum is included in the HS models. The electron temperature, T^, is calculated self-consinstently, and 
thus it is not a model parameter. The parameters marked with 'f ' were kept constant. Fx is the total flux in the unabsorbed model, and Fxo is the model flux 
without the reflection component. The uncertainties correspond to 90 per cent confidence for a parameter, i.e., A^^ = +2.71. 
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allowed range of M. However, thieir constraints can be partially re- 
laxed if the wind is either very hot or cold and only partially ion- 
ized. 

On the other hand, our estimates of M are significantly lower 
than those of van Kerkwijk (1993) and Ogley, Bell Burnell & 
Fender (2001). Van Kerlswijk (1993) have obtained M ^ 4 x 
1O-'']VI0 yr-' based on IR observations of the wind emission, fol- 
lowing the method of Wright & Barlow (1975). Ogley et al. (2001) 
observed IR emission in quiescent and flaring states and found 
M (4-30) X 10-^Mq yr-'. We note that similar discrepancies be- 
tween the M calculated with different methods have been already 
previously observed in both single and binary massive systems 
(IVIoffat & Robert 1994; Nugis, Crowther & Willis 1998; Schild et 
al. 2004). The source of the discrepancy is the usual assumption of 
a smooth wind (Wright & Barlow 1975) when calculating the wind 
IR excess emission. The excess is produced by two-body processes, 
in particular free-free and free-bound, which both the emission and 
absorption are enhanced in a clumpy wind by 1 // with respect to 
the smooth flow, see Section [341 Then, we have found (using the 
formalism of Wright & Barlow 1975 generalized to an arbitrary 
velocity profile) that the net elfect of clumping on measuring the 
mass loss rate using two-body processes and assuming a smooth 
flow is an overestimate by /-'''. This result holds for the measured 
luminosity coming from either the entire wind region (if it is opti- 
cally thin at the observed frequencies) or its outer part of about unit 
optical depth, and it is a generalization of the corresponding result 
of Schild et al. (2004). In the case of a two-phase wind, our mass 
loss rate corresponds to the emission from a smooth flow with the 
rate of [Mj;/(1 - /) + M^/fY'^, see Section below. Given our 
values of / < 10-^, our range of the values of M is consistent with 



van Kerkwijk (1993) and overlap with the lower limit of Ogley et 
al. (2001). 

An additional constraint on M will come from the form of the 
orbital modulation in Cyg X-3 in the HS and SS (shown in H08). 
In this work, we have not attempted such calculations. Some pre- 
liminary modelling of phase-resolved spectra was done in Szostek 
& Zdziarski (2005). 

5.2 Clumpy and inhomogeneous winds 

As discussed above, strong evidence for clumpy winds in WR and 
early type stars comes the discrepancy between estimates of M us- 
ing the IR excess and other methods, which is explained by the en- 
hancement of the two-body emission by /-''^ in a clumpy medium. 
Additonal evidence in single WR stars comes from random vari- 
ability in their flux and polarization (Robert, Moffat & Seggewiss 
1991; Lepine et al. 2000 and references therein). Another piece of 
strong evidence for the inhomogeneities comes from the form of 
the line wings created by frequency redistribution of line photons 
by electron scattering, with clumping reducing the relative contri- 
bution of these wings (Hillier 1991). 

The inhomogeneities can be described as consisting of numer- 
ous dense clumps embedded in a rarefied interclump medium. Ions 
of low and moderate ionization stages are found only in clumps 
while the interclump medium is highly ionized. The cause for for- 
mation of clumps appears to be hydrodynamical instabilities in 
radiatively-driven winds, which lead to shocks, density enhance- 
ments and rarefactions (e.g., Owocki, Castor & Rybicki 1988). 

The evidence for a clumpy wind in Cyg X-3 given in this work 
stems from the form of the X-ray spectra from this object, and it is 
independent of the findings discussed above. A homogeneous wind 
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is found to be strongly ionized and giving rise to an emission excess 
at < 2 keV, not present in the data. The excess disappears when the 
wind is assumed to be clumpy. On the other hand, we still need 
a hot phase to account for the observed Fe edge at ~9 keV. The 
resulting two-phase medium has the relatively low mass fraction in 
the clumpy phase of ~ 1/5-1/3. 

The values of / found by us, < 10"^, are smaller than those of 
other authors. Hillier & Miller (1999) used in their calculation / ~ 
0.1 for WNC (carbon sequence WR) stars. The line profiles fitting 
including electron line wings typically yield / ~ 0.04—0.25 (Morris 
et al. 2000). Schild et al. (2004) found / ^ 0.06 for WR binary 
Velorum. Hamann & Koesterke (1998) obtained / ~ 0.06-0.25. 

We note, however, that those authors assumed empty space 
between the clouds. In our case, the intercloud medium is filled 
with a hot smooth plasma which provides most of the wind mass. 
Then, the density contrast between the two media is > /, and it is 
~ 0.03. Furthermore, we find that for a medium with an arbitrary 
form of inhomogeneities, the enhancement of the wind two-body 
emission is not by /"''^ but by /'"''^, where /' is given by 



/' 



(17) 



In a one-phase medium composed of clouds only, /' = /. However, 
for a two-phase medium, 



r 



(Mh + 



Mll{\-f) + Mllf 



(18) 



which, for our results, yields /' 0.07-0.3, in full agreement with 
the estimates of the clumping factor given above. 

We also point out that the WR star in Cyg X-3, if synchro- 
nized with the binary motion, rotates at a speed comparable to the 
breakup velocity (as noted by Ogley et al. 2001). Namely, its rota- 
tion is 0.12-0.35 of the breakup value for the range of WR models 
considered in Section 1431 This may lead to the wind being com- 
pressed towards the equatorial plane, similarly to the case of Be 
stars (Bjorkman & Cassinelli 1993). This will lead to an anisotropy 
of the wind, not included in our treatment. Furthermore, the wind 
will be focused towards the compact object (Friend & Castor 1982). 
These density enhancement will afi"ect estimates of M from the IR 
emission, and may be partly responsible for the high values ob- 
tained by van Kerkwijk (1993) and Ogley et al. (2001). On the other 
hand, our method of determination of M is sensitive to the density 
along the line of sight, and a wind enhancement away from it will 
be not accounted for in the method. 



5.3 The line spectrum 

In Fig. [131 we plot the 1-12 A (-1-12 keV) SS and HS model 
spectra and identify the strongest emission lines. We plot it in wave- 
length rather than in energy to facilitate a comparison with the re- 
sults of Chandra high resolution spectroscopy (Paerels et al. 2000). 

We find that the lines resolved by Chandra are also the 
strongest lines of our model. However, recombination continua 
which were quite profound in the Chandra observation (in the HS) 
are not visible in the BeppoSAX data. 



5.4 The model of the continuum 

At the best fits, our continua are dominated by Compton reflection. 
In the case of the HS, the the observed shape looks indeed remark- 
ably similar to that of reflection. Recently, a HS Cyg X-3 spectrum 



measured by INTEGRAL was studied by H08, who also found a 
reflection-dominated as the best fit. We discuss below some physi- 
cal implications of this model. 

The reflection component observed in the spectra of both 
XRBs and AGNs may originate in a cold torus surrounding the 
system or an accretion disc. In the case of XRBs, reflection from 
the companion is also possible, though its relative strength is usu- 
ally small. The dominance of reflection means that we do not see 
directly the emitting source, which is obscured, but only its reflec- 
tion from some medium. Such a situation is common in Seyfert 2 
galaxies (e.g.. Matt et al. 1996; Ikebe et al. 2000), where the central 
X-ray source is obscured by a large-scale torus, and we see the re- 
flection from its back side. In XRBs, reflection may take place from 
a torus-like thickened inner accretion flow. A reflection-dominated 
spectrum has also been claimed from another XRB, GX 1+4 (Rea 
et al. 2005). 

Then, Nandra & George (1994) considered reflection from 
cold clouds embedded inside the hot plasma. However, the clouds 
have to be Thomson optically thick for this process to be effective. 
This is not the case in Cyg X-3, where we found that the total op- 
tical depth of the wind along the line of sight is < 1 (Section [4.2t . 
The presence of clumps does not change the average optical depth 
(Osterbrock & Flather 1959), and we can rule out this scenario. 

We have found that the reflecting material is close to neutral, 
and in particular, we see the absorption edge at ~7 keV. On the 
other hand, we found that the Fe K line accompanying reflection is 
fitted at ~6.7 keV, corresponding to highly ionized Fe. This is a se- 
rious problem for the reflection scenario. The observed 7-keV edge 
should be accompanied by a line at 6.4 keV. A line at that energy 
has been observed by ASCA (Kitamoto et al. 1994) and Chandra 
(Paerels et al. 2000). However, its equivalent width was very small, 
e.g., ~60 eV in the latter (HS) data, much smaller than ~1 keV 
predicted in the reflection-dominated scenario. If the observed 6.7 
keV line represented the blue wing of the relativistic disc line with 
the rest energy at 6.4 keV, we should have also seen the K edge 
similarly blueshifted, but we do not. 

We note that a remarkably similar puzzling situation is present 
in some Narrow-Line Seyfert- 1 AGNs, e.g., IH 0707-495 and 
IRAS 13224-3809 (Boiler et al. 2002, 2003; Fabian et al. 2002, 
2004; Tanaka et al. 2004; Gallo et al. 2004). Two classes of physi- 
cal interpretations have been proposed, namely reflection and par- 
tial covering, both by a neutral medium. We consider them here in 
the context of Cyg X-3. 

A reflector in the form of either cold torus or cold wind clumps 
in a natural way produces a 7 keV absorption edge. It also should 
produce neutral iron Ka line but such line is excluded by the data. It 
has been argued that if the reflecting material is in the inner regions 
of the X-ray source and the reflection spectrum inluding the line is 
blurred by Doppler and gravitational redshift effects. The apparent 
edge at ~7 keV is then the blue edge of the relativistically blurred 
iron emission line (Fabian et al. 2004). 

In the other Narrow-Line Seyfert- 1 scenario, the source is par- 
tially covered by a thick cold absorber (Gallo et al. 2004). The par- 
tial absorber does not reduce the soft photon flux as much as a 
complete absorber would do, but, if the absorber is thick enough, 
it can leave its signatures at high energies and produce strong iron 
edge. In addition, similarly as in reflector scenario, if the absorber 
is located close to the compact object, the line will be blurred and 
thus invisible. Such (rather fine-tuned) scenarios appear in principle 
possible for Cyg X-3, but their detailed tests would require further 
studies. 

On the other hand, the wind absorption is very strong in Cyg 
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Figure 13. The best-fit models for the HS (lower, red) and the SS (upper, green) with line identification. The HS model includes the Gaussian line added in 
Section PTTl The SS model also includes a line at =7.6 A (1.6 keV), which we have been unable to identify. The absorption feature in the HS spectrum at =;2.8 
A (4.4 keV) appears to be of instrumental origin. 



X-3 and the form of the spectra corrected for it is rather different 
from that observed, see, e.g.. Fig. [5] Indeed, we find it worrisome 
that the dominance of reflection is found by us in both the HS and 
the SS, whereas the accretion-flow geometry is, in general, depen- 
dent on the accretion rate, with an opticafly-thick disc thickening 
predicted at high rates (e.g., Shakura & Sunyaev 1973; Jaroszynski, 
Abramowicz & Paczynski 1980). We consider it quite possible that 
adding more complexity to our models would change our conclu- 
sion. In particular, a major driver of the enhanced reflection in the 
fit is the presence of the ~ 7 keV edge in the data. That edge may be 
due to partial covering by some weakly ionized additional material, 
e.g., a structure at the outer edge of the accretion disc. 

We point out that the electron temperatures in the Comptoniz- 
ing plasma in the HS obtained by us (Tabled are much lower than 
values typical to the hard state in either black-hole or neutron-star 
binaries. Also, the rather high nonthermal electron fractions found 
by us (Table O are also not typical to the HS. These effects are 
likely to be due to the interaction of the very strong wind with the 
accretion disc, leading e.g. to shocks accelerating electrons out of 
the Maxwellian distribution, see a discussion in H08. 

Thus, the constrains on the shape of the intrinsic continuum 
we have obtained in this paper are relatively limited. We consider 
our present combined wind and continuum study as a first step to- 
wards physical understanding of Cyg X-3. 



5.5 The luminosity 

If we take the results of our spectral fits as a phenomenological ap- 
proximation to the true intrinsic spectra, we can calculate the lumi- 
nosities from the absorption-corrected X-ray bolometric fluxes. Fx 
(Table[3}, assuming isotropy. We also need to correct for the effect 
of the orbital modulation. Our fitted spectra are phase-averaged, 
and the peak fluxes are higher than the average ones by ~ 1/0.75 
(H08). This yields L ^ (0.8^) x lO^s erg s"' (from the HS to SS) at 
d = 9 kpc. This can be compared to the Eddington luminosity for 
a helium system of 2.5 x 10^*(Mx/Mo) erg s"'. For Mx = 20Mo, 
the Eddington ratios are ~2 and ~8 per cent for the HS and SS, re- 
spectively, which are close to those seen in Cyg X-1 (e.g., Zdziarski 
et al. 2002), another high-mass binary wind accretor. This supports 
the black-hole nature of the compact object. We also point out that 
the two other known WR X-ray binaries, IC 10 X-1 and NGC 300 
X-1, which have very similar X-ray luminosities of ~ 10^^ erg s"', 



both appear to contain black holes (Prestwich et al. 2007; Carpano 
et al. 2007a, b). 

If the compact object is a neutron star, these Eddington ra- 
tios become ~25-110 per cent, which are somewhat above those 
seen in the HS and SS of atolls, i.e., weakly-magnetized accret- 
ing neutron-star binaries (Gladstone, Done & Gierlinski 2007). 
However, there are currently no known weakly magnetized neu- 
tron stars in high-mass X-ray binaries. The only hypothetical ex- 
ception is the peculiar system Cir X- 1 ; however, the supergiant na- 
ture of this companion suggested by Jonker, Nelemans & Bassa 
(2007) appears relatively uncertain. On the other hand, the spectra 
of strongly-magnetized neutron-star binaries have an approximate 
bremsstrahlung (with kT ~ 20 keV) form and no detected X-ray 
disc blackbody emission (e.g., Apparao 1994), and thus they look 
completely different from the spectra of Cyg X-3. 

On the other hand, if we take the intrinsic incident spectra im- 
plied by our reflection best-fit models as corresponding to physical 
reality, the assumption of the isotropy has to be dropped as we are 
now looking at the source at a particular inclination at which most 
of the intrinsic flux is not seen. Then, Lx ^ 4ncfFx(i(Sl/4n), where 
Fxo is the observed unabsorbed flux of the direct component (Ta- 
ble [3}. The factor of Cl/4n corresponds to the extreme situation in 
which the actual reflector solid angle is ~ 4n, which assumption 
minimizes the luminosity. This yields Lx > (1.7-16) x 10^^ erg s"'. 
This corresponds to a relatively high Eddington ratio in the HS, im- 
plying, e.g., a high black-hole mass, see a discussion in H08. (Note 
that our incident spectrum in the HS is much harder at £ < 20 keV 
than that fitted by H08 for their reflection-dominated model, which 
results in a higher model luminosity than that found by us.) The 
presence of a neutron-star is unlikely as the Eddington ratio in the 
SS would be then > 5. 



6 CONCLUSIONS 

There are a number of essential features of our model which are 
new in the history of the studies of Cyg X-3. We have used a realis- 
tic model of the WR star, with a heavy small inner core and dense 
wind, and with the abundances typical to helium stars. We have 
used an accurate model of Comptonization together with accurate 
radiative transfer calculations to model the intrinsic X-ray emission 
in the vicinity of the compact object reprocessed by the surrounding 
wind. We have used high quality BeppoSAX data, which are charac- 
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terized by both relatively good spectral resolution and broad-band 
coverage. 

We have used the X-rays emission as a probe of the wind 
properties. We have determined the wind mass-loss rate to be (0.6- 
1.6)x IO^^Mq yr"', which is consistent with mass-loss rates of WR 
stars. We have disproved a previous claim that if Cyg X-3 contains 
a WR star its wind would be completely opaque to X-rays. We have 
also determined the structure of the wind illuminated by the X-rays, 
and found it contains both a hot phase (containing most of its mass) 
and cold dense clumps, occupying less than 1 per cent of the wind 
volume. The hot phase was needed to account for the strong ~9 
keV edge seen in the data, and the cold phase was needed to ab- 
sorb a significant soft-excess emission produced by the hot phase. 
We have also found that the density enhancement in the clumps 
averaged over both phases is consistent with the clumping factor 
found in WR stellar winds using other methods. The lines emitted 
by the wind found in the BeppoSAX data are consistent with the 
high-resolution measurement by Chandra. 

Our fits imply that the observed X-ray continua contain very 
strong Compton reflection components. This would point out to the 
presence of a geometrically thick torus in the accretion flow, which 
would obscure most of the intrinsic X-ray source but allow us to see 
its reflection. The reflecting material is found to be weakly ionized 
(from the presence of an edge at ~7 keV), but no associated fluo- 
rescent Fe Ka (at 6.4 keV) is found. This presents us with a puzzle 
remarkably similar to that encounted in a number of Narrow-Line 
Seyfert 1 galaxies. On the other hand, the finding of the dominance 
of reflection may be an artifact of not including enough details in 
our spectral modelling. In particular, an additional partial covering 
by a weakly ionized medium could be responsible for the ~7-keV 
edge in the data. Based on the bolometric luminosity of the inferred 
intrinsic continua, the compact object in Cyg X-3 most likely a 
black hole. 
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